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Abstract 

We study the infrared singularities associated to 
ultra-soft transverse gluons in the calculation of pho- 
ton production by a quark-gluon plasma. Despite the 
fact that the KLN theorem works in this context and 
provides cancellations of infrared singularities, it does 
not prevent the production rate of low invariant mass 
dileptons to be sensitive to the magnetic mass of glu- 
ons and therefore the rate to be non perturbative. 

LAPTH-756/99, BNL-NT-99/5 



1 Introduction 

It is widely accepted that infrared singularities are 
generally stronger in thermal field theories with 
bosons, compared to their counterparts at zero tem- 
perature. This is due to the singular behavior of the 
Bose-Einstein statistical weight at zero energy, which 



affects massless bosonic field^. As a consequence of 
these stronger singularities, only partial results exist 
concerning their cancellation in the calculation of ob- 
servable quantities in thermal massless theories (see 
01 for instance). So far, there is no general trans- 
lation in the language of thermal field theory of the 
arguments given for this cancellation at T = by 
Kinoshita |Q, and Lee and Nauenberg 

The resummation of hard thermal loops (HTL 
in the following) Q cures partly this problem by 
giving a thermal mass to otherwise massless fields, 
like gauge bosons. Nevertheless, the static magnetic 
(transverse) modes remain massless in this framework 
and may still generate infrared singularities, as ex- 
emplified by the calculation of the fermion damping 
rate |^. In QCD, it is believed that a thermal mass 
for the static transverse modes is generated non per- 
turbatively at the scale g^T, but this mass may be 
too small to be an efficient regulator. 

A particular area where this infrared problem be- 
comes relevant is the thermal production of particles. 
In this paper, we focus mainly on the production 
of photons by a quark gluon plasma. The produc- 
tion rates are calculated as the imaginary part of a 
self-energy diagram evaluated at finite temperature 
1^ , and are expected to be observable quantities that 
should come out finite in a consistent calculation. 

In a recent study iQ, |^, ^, it has been shown that 
2-loop contributions involving the bremsstrahlung 
mechanism overwhelms 1-loop contributions for the 
production of a soft real photon. The insertion of an 



^For a field of mass m, the statistical weight (to be evalu- 
ated on-shell in the real-time formalism) is bounded by T/m. 
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Figure 1: 2-loop dominant contribution to photon 
production. The blob indicates that the gluon propa- 
gator includes the resummation of hard thermal loops. 
The quark propagators include a thermal mass Moo, 
arising from the HTL resummation in the hard limit. 

exchanged gluon in the hard quark loop (see Fig. |l|) 
generates collinear singularities which are power like 
in 2-loop diagrams while they are only logarithmic 
in the 1-loop contributions: as a consequence, when 
these singularities are regularized by the resumma- 
tion of the thermal mass Mao ^ gT on the quark 
propagators, the two-loop diagrams get an enhance- 
ment by powers of g~^ , where g is the strong coupling 
constant. 

The contribution of the diagram of Fig. |^, al- 
though dominated by a soft gluon, is infrared fi- 
nite. In fact, even the contribution of the trans- 
verse gluon is finite in this particular calculation, 
due to kinematical constraints. Indeed, it is trivial 
to see that the two delta functions corresponding to 
the cut quarks 5{P'^ - M^)d{{R + Lf - M^) be- 
come (S(P2 _ M^)S{R^ - M^) in the limit of van- 
ishing L, and that the latter pair of delta functions 
do not have a common support if Mao 7^ 0: for the 
bremsstrahlung process we are considering here the 
energies po and ro have the same sign and hence 
— M^ and — M^ cannot vanish simultane- 
ously, whatever the value of Q^. It is therefore kine- 
matics, via the fermion thermal mass, that prevents 
infrared singularities in this particular topology by 
providing a natural cutoff of order gT on the gluon 
momentum L. This statement was tested in Q by 
studying the limit Moo — *■ 0. A stronger divergence 
was found in the transverse gluon contribution, indi- 
cating that Moo played a role in the regularization of 
this potentially dangerous contribution. 



But this kinematical cut-ofi^ does not apply to ad- 
ditional soft gluons one may insert in the quark loop, 
like in the diagram of Fig. for instance. Indeed, in 




Figure 2: Example of three loop contribution for pho- 
ton production. One of the gluons cannot become 
ultra-soft due to kinematics, and is displayed in bold- 
face. The other, unconstrained, gluon is displayed as 
usual. 

topologies involving more than one exchanged gluon, 
the kinematical argument given above constrains only 
the sum of the momenta of the cut gluons. Therefore, 
we know that L -\- Li cannot vanish, which tells us 
that (for instance) L has a lower bound at the scale 
gT, but Li can still become arbitrary smallQ and this 
leads to an infrared divergence for the cut depicted on 
Fig. H when the ultra-soft gluon is transverse. Indeed 
when compared to the 2-loop diagram, the additional 
gluon provides (i) two coupling constants, (ii) two 
quark propagators, (iii) a set of gluon spectral den- 
sity and statistical weight and (iv) the phase-space of 
the additional gluon. Collecting everything, we can 
estimate by a crude power counting: 

(3 — loop) 

~ (2 - loop) X .92 J d'^Lm^ iLi)p{Li) 
xS{P + Li)S{R + Li) 

~ (2 - loop) X ^ , (1) 

^There is another, symmetric, contribution coming from 
the region of phase space where Li is of order gT and L can 
become arbitrary smaU. 
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where p is the spectral function of the additional 
gluon, where ^ is introduced as a regulator on the 
integral over li . We used the fact that the quarks are 
hard, and mostly on-shellQ because of the cut cross- 
ing the quark loop. It is important to stress here 
that each fermion propagator brings an extra factor li 
in the denominator thus contributing to the infrared 
sensitivity of the above expression. The conclusions 
of this naive power counting are the following: 

• If the additional gluon is longitudinal, its cut- 
off is a thermal mass of order gT , and the 
corresponding contribution is suppressed by one 
power of g compared to the 2-loop one. 

• If the additional gluon is transverse it is natural 
to assume the regulator to be the magnetic mass 
jjL ~ g^T, and we have (3 — loop) ^ (2 — loop). 

Therefore, it seems that if we keep adding trans- 
verse gluons in the quark loop, we generate contri- 
butions that are all of the same order of magnitude. 
This fact is very similar to the argument given by 
Linde for the breakdown of perturbation theory in 
thermal QCD, although in the different context of 
the calculation of the free energy [|lOj . 

There is nevertheless one reason why this power 
counting may be too naive. One should indeed keep 
in mind that this estimate is valid only for a given cut 
through the 3-loop diagram. It does not take into ac- 
count potential compensations that may occur when 
one is summing all the possible cuts. In this paper, 
we are going to study in more detail this possibility, 
and its interplay with a magnetic mass at the scale 

2 Infrared cancellations 

An important feature of the above example is the 
fact that the quark propagators participate in the 
overall infrared divergence of the diagram. In fact, if 
some quark propagators were not becoming singular 
in the IR limit, the diagram would have been finite by 

^ The role played by the small ofT-shellness of the additional 
quarks will be considered later on in this paper. 



power counting. This can be generalized to a topol- 
ogy with an arbitrary number of exchanged gluons 
(but without 3- and 4-gluon vertices)^. Indeed, for 
these topologies, the number of loops L is related to 
the number of gluons Ug by: 

L=\ + ng . (2) 

One of the L loop integrals is an integration d'^P over 
the quark momentum which is hard, and is not con- 
cerned by the IR problem. The remaining L — 1 in- 
tegrals are over the momenta of the Ug soft gluons. 
The fact that L — l = ng tells that even if each gluon 
comes with the singular factor ng{l°)p{L), it is ac- 
companied by a phase space d'^L which is enough to 
make the integral finite. For these topologies, it is 
the quark propagators which are ultimately responsi- 
ble for the IR divergences. Indeed, it is trivial to see 
that if one quark propagator is cut, then some other 
quark propagators become infinite when the gluon 
momenta go to zero. This is the reason why we are 
going to focus on the quarks, and do not care about 
the gluon propagators. 

It is rather easy to show that there are cancella- 
tions at the level of the quark propagators occuring 
in the calculation of the imaginary part of IIJ^. Let 
us illustrate that in the reasonably general situation 
where we have a single quark loop to which are con- 
nected the two external photons, and an arbitrary 
number of internal gluons (but without non-abelian 
couplings, see Fig. ||). To simplify, we detail only the 
lower fermion line, and hide the details of the up- 
per line in a complicated function we do not need to 
specify. 

Since we are going to demonstrate cancellations 
among the various contributions to the dispersive 
part of the photon polarization tensor, it is conve- 
nient to work in the R/A formalism in which 
cutting rules exist that arc both very simple and 
very close to the T — ones Therefore, all 

*In this paper, we are considering only abelian topologies, 
since this is enough for our purpose of studying the interplay 
between the magnetic mass and possible cancellations. Later 
on, we indicate why the arguments given here cannot be ap- 
plied to non-abelian topologies. 

^If the R/A amplitudes are finite, then so do the time- 
ordered ones. Nevertheless, checking the compensations for an 
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Figure 3: Generic configuration for which cancella- 
tions in the infrared sector occur. Again, the gluons 
all include the HTL correction. The gluons do not 
need to he parallel for the argument to work, since 
the gray area hides the details about how the gluons 
are attached to the upper quark line. The index i, 
running from to n, indicates which quark propaga- 
tor is cut on the lower fermion line. 

the contributions to the imaginary part of 11^ can be 
obtained by cuts dividing the diagram into two con- 
nected pieces, each of these parts containing at least 
one external leg. 

By summing over the index of the cut quark on the 
lower fermion line, the contribution of the diagram of 
Fig. H can be written a^ 



A = 



n 

d^p\[{d^Un^{l'l)p{U)TY&ceY{p\p, {Li\) 



i=0 



K,? - Ml 



{P + K, 



Ml 



where the function F hides all the details about the 
denominators on the upper quark line^ (as well as 
the fermionic statistical factors, which do not play 
any role in the following) , and where we have defined 



(4) 



arbitrary number of loops with the rules found in [ p^ would be 
awkward since it is not possible to write all the contributions 
as cut diagrams. If one insists on using^the CTP formalism as 
an exercise, then the rules given in ||l4] are a better starting 
point. 

®In the R/A formalism, we pick the most singular piece for 
each gluon, i.e. (l°)p{Li). Failing to do this, we would get 
a contribution suppressed compared to the 2-loop result. 

^This function does not depend on the position of the cut 
on the lower quark line, but depends on the position of the cut 
on the upper line. 



We can now use the 5 functions to perform explicitly 
the integral over which gives after splitting the 
propagators into positive and negative energy terms 

j d^p[\{d^L.,n,{f,)p{L,) Trace] [] — 



n 



1 



=0 

e = ±l 



(5) 



where we denote uji = ^/{{p + ki)^ + M^) and = 
— fc° + eujj. According to this definition, all the af 
become equal to e-\/(p^ + M^) when the gluon mo- 
menta go to zero. We see that denominators where 
both Qfi's carry the same sign vanish in this limit. 
The IR singularities therefore show up in the vanish- 
ing denominators af — a^. Only the second line in 
Eq. (H) is relevant in the following discussion, and it 
can be compactly rewritten as 



n 

E 

i=0 
{=, = ±1} 



Fial%p,{L.})Y[ 



a. 



(6) 



One can simply observe that for every denominator 
af — a^ with numerator F{af,p, {Li}) appears a de- 
(3) nominator — af with a numerator F{a^ ,p, {Li}) 
(all the other denominators being the same). The 
simple poles therefore cancel trivially. This can be 
extended to the more complicated situation where 
more than two a^s tend to a common value, which 
amounts to prove that these denominators appear in 
a combination that remains finite for any configura- 
tion of the aiS. For that purpose, let us consider an 
expression like^ 



F„. 



n 

E 

4=0 



Fia,)Y[ 



1 



ai — aj 



(7) 



and show that such a quantity is always finite pro- 
vided some regularity property of the function F. 
The shortest way to see that is to notice that F„ 



*At this stage, we can drop all the superscripts since the 
compensations occur in fact for each given set of e^s. 
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is the leading coefficient of the Lagrange polyno- 
mial of degree n that interpolates between the points 
{ai,F{ai)): 



P„(.x)^5]F(a0n- 



(8) 



As such, Fn is finite for every values of the a^'s if 
the function F is n times differentiable. Indeed, if 
several points [ai, F{ai)) collapse into a single point 
of multiplicity m, the Lagrange polynomial has a fi- 
nite limit, and coincides with F and its first m — 1 
derivatives at this point. One can note that this argu- 
ment also applies to insertions of self-energies on the 
quark line. Indeed, the only peculiarity of self-energy 
insertions is that they force several ctiS to have the 
same value even if the Li do not tend to zero; and 
the above proof works for any configuration of the 
aiS, whatever the cause that makes them equal is. 

Although the above argument has been presented 
in the case of a 2-point function, it still works for the 
dispersive part of any n-point function (just attach 
more photons or gluons to the shaded part of Fig. ^). 



transverse gluons which are shielded by a "small" 
magnetic mass. The time-like gluons are shielded 
with a "large" thermal mass of order gT and do not 
lead to any infra-red problems in the context of this 
study. 

We are now going to apply the above considera- 
tions to study explicitly the 3-loop example already 
presented in Fig. |[ As we have already seen in the 
introduction, one of the gluon momenta has a lower 
bound thanks to kinematics, while the other gluon is 
not constrained in the ultra-soft limit. Let us choose 
the gluon on the right (momentum Li) to be ultra- 
soft while momentum L remains soft. This kinemat- 
ical constraint prevents the propagators of the lower 
quarks to become infinite when Li 0: they always 
remain off-shell by some amount controlled by Moa- 
The constraint also simplifies the pattern of cancel- 
lation of infrared divergences in the upper line since 
the propagators of momentum P + L + Li Ri P + L 
when Li — > cannot diverge for the same reason. In 
consequence we need only the two cuts depicted on 
Fig. ^ to get rid of all the zeroes in the denominators^ 
It turns out to be convenient to perform the change 



3 Competition with the mag- 
netic mass 

The result of the previous section is that the prod- 
uct of all the quark propagators remains finite in 
the IR sector when the sum over the cuts has been 
performed. When considering the imaginary part it 
shows that, for a fixed cut on the upper line, there 
is no singularity in the propagators of the lower line. 
The argument should be repeated for the upper line, 
which is made finite by summing over all the ways 
of cutting it. These cancellations between different 
cuts occur within a given topology, and correspond 
to compensation between real and virtual corrections. 
They can therefore be seen as a form of the KLN the- 
orem. 

We point out some differences with the usual ver- 
sion of the KLN theorem at zero temperature where 
the exchanged gluons are bare gluons. At finite tem- 
perature we deal with resummed gluons and the most 
dangerous divergences arise when cutting space-like 



,- a p ,b 




Figure 4: Notations used for the 3-loop contributions. 



of variable P + Li ^ P in the second contribution, 
as indicated on Fig. |^. Moreover, we neglect Li in 
front of L whenever these two impulsions appear to- 
gether (like in Li + L fi: L). With these notations 
and approximations, one can readily see that the ad- 
ditional ultra-soft gluon brings the following factors, 
to be multiplied by the 2-loop integrand (i.e. the in- 
tegrand for the diagram of Fig. |l|), respectively for 

^We also have to take into account the other possibility with 
Ll soft and L ultra-soft which is easily done by multiplying the 
final result by a factor 2. 
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the cuts (a) and (&): 



Fa,b = r 



(2P,)(2i?,)P;''(ii) 



((P±Li)^-M^)((i?±ii) 



,(9) 



where P^'' is the transverse projector, and is the 
transverse spectral function of the gluon (only the 
transverse mode of the ultra-soft gluon is relevant, 
the other gluon can be transverse or longitudinal, 
as in the 2-loop calculation). The factors of form 
Pa/ {{P ± Li)^ — M^) come from the fcrmion prop- 
agators in the hard momentum approximation. We 
can further simplify these factors by evaluating them 
at P^ = since they are to be multiplied by the 2- 
loop integrand which contains a S{P^ — M^). Notic- 
ing that (see eq. (46) in [|j) 

PaRpP^PiL^) « -pr{\ - cos^ Oi) , (10) 

where 9i is the angle between the vectors p and Zi, 
it is easy to perform the angular integral explicitly, 
which gives: 



1 - cos^ 6i 





r A 1 




( 




In 




L 4rZi J 





± 



{{P±L,f 

TT 

prll 
7r(l - 
phA 

-In 



Ml){{R±Li[ 
2rli{x + 



2 _ ]\,p 
1)± A 



2rli{x- 1)± A 



In 


x + 1 


x-1 





2rli(x + l)±A 



(11) 



2rh{x- 1) ± A 

where we denote x = l^h and A = 2P-Q + Q^. The 
general result established in the previous section says 
that the sum over the cuts for this quantity should 
have a finite limit when ^ at fixed x (i.e. when 
the 4 components of Li tend to zero). One can see 
that this is not the case on the contributions of the 
individual cuts, since we have: 



1 - cos2 01 



{{P±L,r-Ml)iiR±L^r-Ml 



± 



7r(l 



phA 



In 



x+1 



± 



Anx 
PhA 



(12) 



(where we have neglected terms regular in li) but oc- 
curs trivially for the sum of the two cuts since the sin- 
gular (l^^) part drops out. The interesting quantity 
is in fact the scale at which this cancellation occurs. 
The only dimensional quantity to which li can be 
compared with in the above expressions is ^ A/r. 
Carrying out the integration on x, li we find the fol- 
lowing finite result: 



Fa. 



+b 



g'^Tr 
2A 



(13) 



When we add a magnetic mass /i into the game, 
we have two regularization mechanisms in competi- 
tion, and we must compare the scales at which they 
operate, the most efficient being the one that has the 
largest scale. A convenient way to introduce^ the 
magnetic mass is via the following sum rule 



+ 00 



dx p^{li,x) 
27r X 



B{x) 



B{0) 



(14) 



which is a reasonable approximation if B is not singu- 
lar at a; = and does not increase too much when x 
becomes large. Adding the two cuts and approximat- 
ing Tig {li) « T/hx, we find that the 2-loop integrand 
is to be multiplied by the factor 



Fa. 



+b 



+ 00 

g^Tr f hdl 

TT^A 



In 



A + 2rli 



A - 2rh 



A 




A 


+ 2rli 










- 1 






A 


-2rh 



} • (15) 



We can now give analytical limits for this factor in 
two cases. If the magnetic mass dominates over the 
scale II defined above, we have 



Fa 



+b 



_ 9 



1 



(16) 



while in the opposite limit {l^ <^ ii) we recover 
Eq. (©. 

^''This does not tell much about the way the magnetic mass 
enters in the gluon propagator: it just tells that the self-energy 
of the transverse gluon satisfies Tl^{x = 0) = fi^, which is a 
definition analogous to that of the Debye screening mass. 
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Let us mention that if the additional gluon in Fig. Q 
is longitudinal a similar calculation as above goes 
through where fi is to be replaced by the Debye mass 
of order gT, so that the factor Fa+b is of order .g ^ 1. 
The conclusion is therefore the following: 

• If the magnetic mass is the relevant regula- 
tor, then the 3-loop diagram gives a "correction" 
of order one to the prefactor of the 2-loop re- 
sult, and this is likely to be the case for higher 
loop corrections also. In this regime, the pho- 
ton production rate is sensitive to the magnetic 
mass, and one must resum an infinite series of 
diagrams. 

• If on the contrary, the scale A/r dominates or 
if the regulator on the gluon propagator is the 
Debye mass, then the 3-loop diagrams lead to a 
negligeable correction to the 2-loop result, and 
it is expected that this hierarchy will be valid 
between 4-loop and 3-loop,... 1^. 

It remains now to make a bit more explicit the 
comparison between A/r and /i. For that purpose, 
let us recall that the expression of A = 2P ■ Q + 
can be taken from the 2-loop calculation [Q: 



'pqo 



1 — cos 9 



Ml 



2p' 



(18) 



where 9 is the angle between p and q, and where^ 

(19) 



—pr 



Due to the extremely singular nature of the integral 
over 6* 1^, |j, the order of magnitude of A is qoM^g/p. 
Taking into account the fact that r ~ p + qq and 
p ~ T, we divide the {q,qo) plane in two regions 
where respectively fj, or A/r dominates (see Fig. ^). 
In region I, i.e. roughly the region of small virtual- 
ity, the magnetic mass is the most important regu- 
lator and, as seen above, the 3-loop contribution is 
as large as the 2-loop contribution: the production 
mechanism becomes non-perturbative and resumma- 
tions should be considered. Increasing the virtuality 

^^The formula for has been extended here to hold for 

hard slightly virtual photons |l5| . 



Il-b 



Il-a 



_g2T 



Figure 5: Comparison of fi and A/r in the {q,qo) 
plane. In region I, the magnetic mass is the most im- 
portant regulator, the 3-loop contribution is as large 
as the 2-loop one. In regions Il-a and Il-h, the 
magnetic mass is a suhdominant regulator, and the 
higher-loop contributions are suhdominant. In region 
Il-h, L can become hard and some of our approxima- 
tions become invalid (see J^). 

of the photon increases A (see Eq. ([l9|)) which 
eventually becomes the dominant cut-off and one en- 
ters region II where the infrared sensitivity of 3-loop 
(and higher loop diagrams) becomes suhdominant. 

In region Il-b, the production mechanism receives 
a contribution from large gluon momentum and the 
approximations done in the previous calculations may 
become incorrect. However, it was found by an ex- 
plicit calculation that the expressions valid in re- 
gion Il-a could be safely extrapolated to the case of 
soft virtual photon at rest {q — 0) in region Il-b. 
However, there is the possibility that 3-loop diagrams 
give important contributions in region III due to hard 
gluon exchanges, but this is a different story. 

4 Summary and outlook 

In this paper, we have considered higher order correc- 
tions to photon production by a quark-gluon plasma, 
with emphasis on the infrared singularities due to 
ultra-soft transverse giuons. When one sums over all 
the possible cuts through a given abelian topology, 
there are cancellations that prevent the quark prop- 
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agators from becoming infinite, making the higher- 
loop corrections infrared finite. 

The generalization of this result to non-abelian 
topologies is not straightforward. Indeed, the above 
argument is flawed already at the stage of counting 
the number of gluon propagators. Loosely speaking, 
because of the 3- and 4-gluon couplings, the number 
of gluons can be larger than the number of indepen- 
dent momenta. The loop counting for QCD gives for 
the number of gluon integrations 

L - I = Ug - + Hi) < Hg , (20) 

where and are respectively the number of 3- 
and 4-gluon vertices. As a consequence, unlike in the 
abelian topologies for which L — 1 = Ug, the gluon 
propagators play an active role in the KLN compen- 
sations (if any), and much more elaborate arguments 
are required. As a side note, it is also impossible to 
apply in the thermal case the heuristic power count- 
ing argument of Poggio and Quinn since some 
gluons have an infrared count of —3 (instead of —2 
at T = 0) because of the Bose-Einstein statistical 
weight. 

We have also considered the possible competition 
of this mechanism of regularization with a hypothet- 
ical magnetic mass of order g^T. The result is that 
there exists a region in the (g, go) plane (roughly 
speaking, for photons of small invariant mass) where 
the magnetic mass is the dominant regulator, and 
where the higher-loop corrections are as large as the 
two-loop result, and should be resummed. 

Further extensions of this work include the treat- 
ment of non-abelian topologies, the study of resum- 
mations in the region where the magnetic mass domi- 
nates, the extension of the 2-loop calculation for hard 
massive photons, and also the study of how the finite 
lifetime of quarks in a plasma affects the coUinear 
singularities found at 2-loop. 
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